Is it possible to observe experimentally a metal-insulator transition in ultra cold 

atoms? 
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It has been recently reported [TrI ] that kicked rotators with certain non-analytic potentials avoid 
dynamical localization and undergo a metal-insulator transition. We show that typical properties 
of this transition are still present as the non-analyticity is progressively smoothed out provided 
that the smoothing is less than a certain limiting value. We have identified a smoothing dependent 
time scale such that full dynamical localization is absent and the quantum momentum distribution 
develops power-law tails with anomalous decay exponents as in the case of a conductor at the metal- 
insulator transition. We discuss under what conditions these findings may be verified experimentally 
by using ultra cold atoms techniques. It is found that ultra-cold atoms can indeed be utilized for 
the experimental investigation of the metal-insulator transition. 
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PACS numbers: 72.15.Rn, 71.30.+h, 05.45.Df, 05.40.-a 

The study of a quantum particle in a random poten- 
tial Q is one of the cornerstones of modern condensed 
matter physics. In its simplest form, namely, a free spin- 
less particle in a short-range disordered potential with 
no interactions at zero temperature, the combination of 
the one parameter scaling theory , the supersymmetry 
method 31 and numerical simulations has led to the 
following picture: In two and lower dimensions destruc- 
tive interference caused by backscattering produces expo- 
nential localization of the eigenstates in real space for any 
amount of disorder. As a consequence, quantum trans- 
port is suppressed, the spectrum is uncorrelated (Pois- 
son) and the system becomes an insulator. In more than 
two dimensions there exists a metal insulator transition 
(usually referred to as Anderson transition (AT))for a 
critical amount of disorder. By critical disorder we mean 
a disorder such that, if increased, all the eigenstates be- 
come exponentially localized. For a disorder strength be- 
low the critical one, the system has a mobility edge at a 
certain energy which separates localized from delocalized 
states. Its position moves away from the band center as 
the disorder is decreased. Delocalized eigenstates, typical 
of a metal, are extended through the sample and the level 
statistics agree with the random matrix prediction for the 
appropriate symmetry. In three and higher dimensions 
the AT takes place in a region of strong disorder only ac- 
cessible to numerical 0, U simulations. Typical features 
of the AT include: 

1. The spectrum of the Hamiltonian is scale invari- 
ant |(J, namely, any spectral correlator utilized to de- 
scribe the spectral properties of the disordered Hamil- 
tonian does not depend on the system size. The spec- 
tral correlations at the AT, usually referred to as critical 
statistics 0, Q , are intermediate between that of a metal 
and that of an insulator. 



2. Anomalous scaling of the eigenfunction moments, 
V q = f d d r\ip(r)\ 2q cx L~ D " { - q ' 1) with respect to the sam- 
ple size L, where D q is a set of exponents describing the 
AT. Eigenfunctions with such a nontrivial (multi) scaling 
are usually dubbed multifractals [f| (for a review see 

3. Quantum diffusion is anomalous at the AT. In 
the metallic limit, up to small weak localization correc- 
tions, the density of probability is Gaussian-like and the 
dynamics is well described by a Brownian motion. How- 
ever, as disorder increases, localization effects become 
important and quantum diffusion slows down. The den- 
sity of probability develops power-law tails with a decay 
exponent depending on the spectrum of multifractal di- 
mensions 0. 

Unfortunately the experimental verification of the AT 
is a challenging task. In the context of electronic sys- 
tems is extremely hard to disentangle effects caused by 
short decoherence times, electron-electron interactions 
and phonon-electron interactions from destructive quan- 
tum interference and symmetry, supposed to be the main 
ingredients driving the AT. 

In recent years ultracold atoms in optical lattices 01 
has been utilized to model certain solid state physics sys- 
tems. Generically, in these experiments a very dilute 
almost free gas of atoms (Cs and Rb) is cooled up to 
temperatures of the order of tens \xK and then inter- 
acts with an optical lattice. In its simplest form, the 
optical lattice consists of two laser beams prepared in 
such a way that the resulting interference pattern is a 
stationary plane wave in space. The laser frequency is 
tuned close to a resonance of the atomic system in or- 
der to enhance the atom-laser coupling but not too close 
to avoid spontaneous emission. In this limit the system 
laser-atom can be considered as a point particle in a sine 
potential, namely, the quantum pendulum. Additionally 
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if the laser is turned on only in a series of short periodic 
pulses the resulting system is very well approximated by 
the so called quantum kicked rotor (see [ll| for a review) 
extensively studied in the context of quantum chaos, 

H=p 2 /2 + kcos(q)J2 S ( t - Tn )- (!) 

n 

The classical motion of this system is diffusive in momen- 
tum space. For short time scales, quantum and classical 
motion agrees. However quantum diffusion is eventually 
suppressed due to interference effects and eigenstates are 
exponentially localized in momentum space. This coun- 
terintuitive feature, usually referred to as dynamical lo- 
calization 01 j was fully understood 0] after mapping 
the kicked rotator problem onto an short range ID disor- 
dered system where localization is well established. The 
first direct experimental realization of the kicked rotor 
was reported in Ref. 0|. As was expected, the output 
of the experiment (the distribution of the atom momen- 
tum and the energy diffusion as a function of time) fully 
agrees with the theoretical prediction of dynamical local- 
ization |l3j . Finally we remark that, after the pioneering 
work of Ref.fiol]. many other aspects of the physics of 
a quantum kicked rotor as the effect of noise and dissi- 
pation have also been investigated |14| by using similar 
experimental settings. 

The above results do not depend on the exact details of 
the potential but only on its ability to produce classical 
chaotic motion. The situation is different if the potential 
is not smooth. Recently 0] it has been reported that 
a kicked rotor could avoid full dynamical localization if 
the smooth sinusoidal optical potential is replaced with 
a generic potential with a logarithmic or step like singu- 
larity. It was found that, for these potentials, the kicked 
particle has striking similarities with a free particle in a 
disordered potential at the AT. Thus level statistics are 
given by critical statistics, eigenfunctions are multifractal 
and quantum diffusion becomes anomalous. 

A natural question to ask is whether this non- 
analytical kicked rotor can be realized in experiments. 
If so, this would be an ideal setting to test the physics 
of the AT. Obviously, in experiments the singularity can 
only be approximated. For instance, an optical lattice 
potential with an approximate step-like singularity can 
be produced either by a holographic mask with 
precision a or by adding a limited number of Fourier com- 
ponents. In both cases the potential is smooth on suffi- 
ciently small scales ~ a . That means that for momenta 
Pd ^ h/cr and times td sufficiently long, the microscopic 
smoothness of the potential is at work and standard dy- 
namical localization should be observed. On the other 
hand, for momenta p c and times t c sufficiently short, 
classical and quantum results should coincide. In be- 
tween these two scales typical properties of the AT are 
observed. The aim of this paper is twofold on the one 
hand we seek to determine in what window of a the AT 
is observed. On the other hand we examine whether this 
range is already experimentally accessible by using ultra 



cold atoms in optical lattices. 

The organization of the paper is as follows: In the 
next section we introduce a kicked rotor with two dif- 
ferent smoothed out versions of a step potential. Then 
we evaluate the rate of energy diffusion and the full mo- 
mentum distribution. Finally we establish the minimum 
smoothing required to observe the AT and whether an ex- 
perimental verification is realistic with the current, state 
of the art, ultra cold atom techniques. 

I. THE MODEL AND OBSERVABLES 

We investigate a kicked rotor in 1 + 1D with a smoothed 
step-like potential, 

H = ?j + V 1 , 2 ( q )J2Ht-nT) (2) 

n 

with q € [— 7r, 7r). We consider the following two poten- 
tials, 

Vi(q) = Si ((£ + q)/a) + St ((| - q)/a) (3) 

where Si(q) = Sln ^ dt is the sine integral function, 
and 

M 

V 2 {q) = ^ f{m)cos(mq) (4) 

m— 

where f(m) is the discrete Fourier transform of the bare 
step like potential V(q) = tt for \q\ < ir/2 and zero oth- 
erwise. In both cases for a — > (a = 1/M in the lat- 
ter case) we recover the bare step-like potential investi- 
gated in |Tfij . Obviously there are infinitely many ways 
to smooth a singularity, we have chosen the above two 
due to similarities with the experimental situation. Thus 
V\{q) represents an optical lattice with square- wave in- 
tensity profile as produced by an array of fine slits or a 
holographic mask [l^. The other potential Vz{q) pro- 
duces an approximated step-like shape by adding a lim- 
ited number of Fourier components. We remark that re- 
sults for Vi(q) and ^2(9) are hardly distinguishable, both 
are smooth and oscillatory on scales of the order of a. 
Numerically it is a little easier to simulate V\{q) so we 
will stick to it for the numerical calculations. 

We analyze both the classical and the quantum motion 
of the above Hamiltonian. The classical evolution over 
a period T is dictated by the map: p n +i = Pn — 9 ^q"^ > 
q n +i = In + Tp n+ \ (mod 2ir). By smoothing the step 
potential the classical force has a well defined classical 
limit for any finite a. 

The quantum dynamics is governed by the quantum 
evolution operator hi over a period T. Thus, after a pe- 
riod T, an initial state ipa evolves to ij)(T) = Uipo = 

-ip 2 T iV(q) 

e 2h e h tpo where p and q stands for the usual mo- 
mentum and position operator. Our aim is to evolve a 
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FIG. 1: (Color online) Comparison of quantum and classical 
energy diffusion versus time for different smoothings: a — 
5 x 10~ 2 (a), a = 10" 2 (b) and a = 10~ 3 (c). Quantum initial 
conditions, IV'(O)) = |0), were chosen to mimic its classical 
counterpart. For t < t c quantum and classical diffusion rates 
are similar. For td ~ 2/er the quantum energy diffusion gets 
saturated due to destructive interference. In between these 
two scales the system behaves as a disordered conductor at 
the AT. 



given initial state to a certain time nT . This is equiva- 
lent to solving the eigenvalue problem U^ n = e~ tKn / h 'fy n 
where ^ n is an eigenstate of U with quasi-eigenvalue K n . 
In order to proceed we can express the evolution opera- 
tor (m\U\n) — U nm in the basis of momentum eigenstates 



Umn — 



2tt 



dqe 



iq(m—n) — iV(q)/h 



(5) 



We remark that in this representation, referred to as 
'cylinder representation', the resulting matrix U nm is uni- 
tary exclusively in the N — > oo limit. This is certainly a 
disadvantage since besides typical finite size effects one 
has also to face truncation effects, namely, the integral 
of the density of probability is not exactly the unity and 
eigenvalues are not pure phases (e~ ie ") as expected in 
a Unitary matrix. Moreover the diagonalization of a 
generic non Unitary matrix is numerically much more 
demanding. 

These difficulties can be circumvented by changing rep- 
resentations in each quantum iteration step, a technique 
extensively adopted in quantum kicked rotator studies. 
First, we express a given state ip in position representa- 

iV(q) 

tion, so that it is straightforward to get tp — e h 4>, 
the state just after the kick. Next, we express ip' in 
the angular momentum representation by using the fast 
Fourier transformation (FFT) algorithm to facilitate the 

-ip 2 T 

calculation of e ^ tp . Since no matrix diagonalization 
is involved in this scheme, the computation is quite fast 
and the effective dimension of the state vectors is as large 
as 10 s . As a result the truncation effects mentioned above 
can be safely neglected. We recall that this method al- 
lows to resolve the potential with a precision of 10~ 8 , 
four order less than the minimum cr(10 -4 ) investigated. 
Such degree of precision is a necessary requirement to de- 
termine the effect of a small a in the quantum transport 
properties of the model studied. 

Analytical results for the above model can in princi- 
ple be obtained by mapping Eq.J2J onto an ID Anderson 
model. This method was introduced in [13j for the case 
of a kicked rotor with a smooth potential. We do not re- 
peat here the details of the calculation but just state how 
the ID Anderson model is modified by the non-analytical 
potential. It turns out that the classical non-analyticity 
induces long-range disorder in the associated ID Ander- 
son model. If the kick strength is sufficiently large the di- 
agonal part of the Anderson model is pseudo-random and 
the off-diagonal one decays as U r ~ 1/r with r the dis- 
tance from the diagonal. This Anderson model is similar 
to the one studied in [5| which is solved by using the su- 
persymmetry method. In general, according to Ref.[{|, a 
1/r decay in ID is a signature of an AT. For the potential 
V\ } 2 above, it is straightforward to show that U r ~ 1/r for 
r ^ I /a and U r ~ e~ ar for r > l/cr. Consequently we 
expect to observe AT like behavior for small momentum 
and then eventually recover the results of the sinusoidal 
potential, namely, exponential localization in momentum 
space. For further details of the analytical approach we 
refer to |l5| . 

We are mainly interested in observables related to 
transport properties as the density of probability and the 
rate of diffusion. 
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The density of probability (both classical P(p,t) and 
quantum P q (p,t)) of finding a particle with momentum 
p after a time t for a given initial state \ip(0)) = |0). 
P q (p,t) = Pg(k,t) = \(k\(j>(t)}\ 2 with p = kh. In all 
calculations we set % = 1. The classical P(p,t) is ob- 
tained by evolving the classical equation of motion for 
2 x 10 7 different initial conditions with zero momentum 
p = and uniformly distributed position along the in- 
terval (— 7r,7r). We would like to emphasize our results 
do not depend on the initial conditions. For instance, 
we have checked that similar results are obtained if the 
initial conditions of Ref.[l6l| are considered. 

We also examine the second moment of the probabil- 
ity distribution, namely, the energy diffusion (p 2 (t)) = 
J Q dpp 2 P(p, t) as a function of time. 

We recall our aim is to find out whether the transport 
properties are compatibles with those of a disordered con- 
ductor at the AT and how they are affected by the short 
distance differentiability of the potential. For the sake 
of completeness let us briefly summarize the predictions 
for both a kicked particle in a smooth potential and a 
disordered conductor at the AT. 

For a kicked rotator with a smooth potential, it is well 
established that initially (up to a certain time t c ) both 
classical and quantum probabilities are Gaussian like and 
the diffusion in momentum is normal, namely, a standard 
Brownian motion. For longer times the classical density 
of probability is still that of a normal diffusion process. 
However P q (p, t) become exponentially localized and en- 
ergy diffusion stops (p 2 (t)) ~ cons. These are typical 
signatures of dynamical localization. 

At the AT, up to a certain t c , agreement is also ex- 
pected between the classical and quantum predictions. In 
the case of a disordered conductor the classical dynamics 
is obviously well described by a Brownian motion. How- 
ever, for t > t c , the diffusion becomes anomalous, the 
quantum density of probability develop power-law tails 
in space (localization in a disordered conductor occurs 
in real space) and time with exponents related to the 
multifractal dimensions of the eigenstates . The rate 
of diffusion is in some cases still similar to the one cor- 
responding to normal diffusion (p 2 (t)) = D quan t though 
the quantum diffusion constant D quan is typically lower. 
This suggests that, at the AT, destructive interference 
is still at work but it is not sufficient to fully localize 
the particle. In our case we also expect agreement be- 
tween classical and quantum results up to a certain time 
i c (cr). Additionally, since the potential is differentiable 
for distances smaller than the smoothing a, we expect 
that there exists a td{o~) such that for i > standard 
dynamical localization becomes dominant. 

Typical features of the AT transition are thus observed 
in our model only if t c <C td- It is unclear for what range 
of rx, td S> t c and whether these values of a can be reached 
experimentally. We answer these questions in the next 
section. 
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FIG. 2: (Color online) Classical density of probability distri- 
bution P(p',t) for a = I0~ 3 and p = p/V2Dt. (a) Region of 
anomalous diffusion, t < 10, (b) Crossover from anomalous to 
normal diffusion 10 < t < td ~ 2/ a and (c) Normal diffusion 
t > td- The results were obtained after averaging over 2 x I0 7 
initial conditions. 



II. RESULTS 

For the sake of clearness we first enunciate our main 
conclusions: 

1. For a < 0.05 we have observed typical signatures of 
an AT in a broad region of times t c > i > 

2. The quantum-classical breaking time t c decreases 
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weakly with a. In the range of er investigated it never 
goes beyond a few kicks. By contrast, the time scale 
signaling the beginning of full dynamical localization, due 
to the differentiability of the potential, increases as a 
decreases, td~2/a. 

3. We argue that the above range of parameters is 
accessible to experimental verification. By using holo- 
graphic mask techniques one can reach up to er ~ 0.01 
|17|. On the other hand coherence in ultra cold atoms is 
maintained well beyond 1000 kicks. Consequently the AT 
can be investigated by using ultracold atoms in optical 
lattices. 

We have computed (see details in previous section) 
the quantum and classical density of probability for the 
Hamiltonian Eq. [21 with potentials given by Eq|3] and a 
variety of smoothings a € [0.1, 10~ 4 ]. 

Our first task is to determine t c and td as a function 
of er. These time scales can in principle be calculated by 
using different observables. Qualitatively all observables 
should provide the same physical picture. However the 
numerical value of t c and td may depend on the observ- 
able considered. For the sake of simplicity we estimate 
these time scales by looking at the the rate of energy 
diffusion (p 2 (t)). 



A. Energy diffusion 

As was mentioned previously we used as initial condi- 
tions p — and random positions. Obviously only initial 
conditions in the narrow region [— 7r+cr, it— a] get a sizable 
kick. Thus even after several kicks there is a high proba- 
bility that the system stays in the region p = 0. In order 
to show that our results are not sensitive to initial con- 
ditions and stable under perturbations we have added a 
weak noise V(q) = fcsin(g) with k = 1. We have checked 
our results do not depend on k provided k <C 1/cr. Ob- 
viously for k ~ I/cr the effect of the pseudo-singularity 
is obscured by the noise strength. In the classical case 
(see Fig.l ) (p 2 (t)) increases linearly with time. The de- 
pendence of the diffusion coefficient on a is well approxi- 
mated by 0.5/fi. This is consistent with the the analytical 
prediction resulting from the random phase approxima- 
tion . In the quantum case (see Fig. I ) we distinguish 
three different regions. In a first stage (t < t c < 10) 
the quantum averaged energy moves around its classical 
counterpart. t c depends weekly on a, it decreases as a 
does. For longer times t c < t < td, the diffusion is still 
linear, (p 2 (t)) « D quan t, with D quan ~ 0.2/cr. Although 
it has the same dependence on tr, it is smaller than in the 
classical case. This suggests that quantum interference 
effects slow down the classical diffusion. A similar feature 
has been found in a disordered conductor at the AT 0. 
This stage lasts up to td ~ 2/ a. For longer times stan- 
dard dynamical localization due to the differentiability of 
the potential takes over and diffusion stops. 

We recall that linear energy diffusion is only a neces- 
sary condition for normal diffusion. In general the infor- 
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FIG. 3: (Color online) Quantum density of probability distri- 
bution P q (p',t), p' = p/V2Dt and a = 1CT 3 for three differ- 
ent regions: (a) t < t c < 10, agreement between classical and 
quantum results, (b) t c < t < td, typical properties of an AT 
are observed, and (c) t > td, standard dynamical localization 
for |p| > y/2Dtd- In the insets we present the same results 
in a linear-log scale. P q (p',t) ~ 1/p 2 for t < t c and \p\' < 1, 
and P q (p',t) ~ l/p a with a = 1.1 ± 0.2 for t > t c . As initial 
condition we used |i/>(0)) = |0). 



mation obtained from the knowledge of a few moments 
of the distribution is not sufficient to fully characterize 
the classical motion. Thus the second moment may be 
(p 2 (t)) ~ t but this by no means assures that the density 
of probability is Gaussian- like |20| as for normal diffusion. 
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We show below that this is the case in our model. 



B. Density of probability 

In the classical case we distinguish two different 
regimes separated by a broad crossover region (see den- 
sity of probability in Fig. 2): First, for short time scales 
(a few kicks) and \p\ < c(a)y/2Dt the diffusion is anoma- 
lous. P(p, t) ~ p~ a with a ~ 2 and c(a) ~ 1 slightly 
increases as a is decreased. For such a short time scale 
the classical system does not feel the differentiability of 
the potential. 

For longer times but t < td, we observe a gradual 
crossover from anomalous to normal diffusion. For small 
momentum the density is still non Gaussian as the effect 
of the pseudo non differentiability is still important. As 
time approaches td, the central (small momentum) non- 
Gaussian region becomes smaller and smaller. Mean- 
while, the outskirts bend down and a Gaussian-like be- 
havior typical of normal diffusion is observed. Finally, 
for t > td, P(p,t) is well approximated by a Gaussian 
distribution. These regions have been observed for all a 
of interest. 

In the quantum case three regimes are distinguished 
(see Fig. 3): 

1. t < t c and \p\ < c(a)V2Dt with D w 0.5/er (c(cr) 
increases slightly as a decreases). The classical and the 
quantum probability agree in this region. The scale t c 
depends weakly on a; it decreases as a does. We recall 
that our system has not a well defined classical-quantum 
correspondence in the limit a — > 0. In both cases the 
diffusion is anomalous P{p,t) — P q (p,t) ~ 1/p 2 - We re- 
mark P q (p, t) is calculated by summing the probability of 
falling in a bin of width Ap' — 0.03 (so does for classical 
P(p,t)), hence it is in fact a coarse-grained result where 
part of quantum fluct uatio ns have been suppressed. 

2. For \p\ < c(a)V2Di but t d > t > t c . The quan- 
tum probability P q {p,t) ~ l/p a develop an power-law 
tail with an exponent a < 2 (see Fig. 3) typical of anoma- 
lous diffusion. The exponent a does not depend on er, in 
all cases we have found a ~ 1.1 ± 0.2. This is a clear 
signature of an AT. We remark that, in agreement with 
previous results from the energy diffusion, the quantum 
decay is slower than the classical one. Quantum interfer- 
ence slows down the motion but it is not enough to fully 
localize the particle. 

3. For \p\ > c(a)y/2Dt and t > t d , P q (p,t) decays 
exponentially. This is an indication of full dynamical 
localization due to the differentiability of the potential. 

From the above we can affirm that in order to observe 
typical features of an AT in the transport properties of 
our system, t c and td must be well separated, namely, 
td 3> t ^> t c . As is shown in Fig. 3 and Table I, this 
occurs provided that a < 0.05. Thus for an experimental 
verification of the AT in cold atoms one has to manage to 
produce a bare step-potential up to corrections of order 
a < 0.05 
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TABLE I: Time scales t c and td for various values of a 



C. Experimental verification 

A natural question to ask is what is the minimum val- 
ues of a that it can be reached in experiments. Specifi- 
cally, we wish to determine, for instance, the maximum 
number of terms in V\ ( q) that can be included experimen- 
tally. In principle [l7| it is an challenging experimental 
task to realize optical potentials with high slopes involv- 
ing higher optical harmonics of the laser beam. The prob- 
lem is that, for instance for Cs, the fourth harmonic is 
already in the vacuum UV, and difficult to produce. Ad- 
ditionally higher order harmonics are not resonant with 
the atom and need a much stronger intensity. Thus it 
seems extremely hard to go beyond the first few harmon- 
ics. Another option could be to use a kicked rotor with 
a smooth potential and three incommensurate frequen- 
cies. According to the results of [2l]], this model can be 
mapped it onto a 3D Anderson model which is supposed 
to undergo an AT for a specific value of the coupling con- 
stant. However in more than one dimension there is no 
clear evidence that this mapping is really accurate. For 
instance, the critical exponents at the AT are very dif- 
ferent from the one found in the kicked rotor with three 
incommensurate frequencies |22j . 

A more promising alternative is to use a holographic 
mask to give a square- wave intensity profile |18j. This 
technique combined with the recent introduction of spa- 
tial light modulators permit the production of a very 
broad range of intensity pattern which act as a effective 
spatial potential for atoms. Unlike the previous method 
the sharpness of the edges would be limited by diffraction 
effects of the order of the wavelength With the cur- 
rent techniques the potential of Eq|3| could be produced 
in a window a > 10~ 2 . On the other hand quantum co- 
herence in cold atoms is lost after a few thousands kicks. 
The experimental bounds are thus within (see above) the 
theoretical limits and, as a consequence, the AT can be 
studied by using ultra cold atoms in optical lattices. 

In conclusion, we have explicitly shown that kick ro- 
tors with a singular but slightly smoothed potential still 
have similar transport properties that those of a disor- 
dered conductor at the AT provided that the degree of 
smoothing is weak enough. The utilization of ultra cold 
atoms in optical lattices offers the opportunity to inves- 
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tigate Anderson localization in general and the AT in 
particular in a setting free from many of the inconve- 
niences that have plagued other experimental studies of 
the AT in the context of condensed matter physics. 
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